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Abstract 

The  model  used  Is  that  of  two  fluids  of  Infinite  deptn,  with 
the  Interface  Initially  In  the  form  of  a  sine  wave  with  amplitude 
small  compared  to  wave  length.  The  fluids  are  considered  Incom¬ 
pressible,  and  only  the  linear  terms  In  the  equations  of  hydro¬ 
dynamics  are  used.  The  first  four  sections  discuss  the  effects 
of  surface  tension  and  viscosity.  The  fifth  gives  a  few  numerical 
results  to  Illustrate  the  main  points  of  the  preceding  sections.  ( 


EFFECTS  OF  SURFACE  TENSION  AND  VISCOSITY  ON 
TAYLOR  INSTABILITY 


Richard  Bellman 
Ralph  H.  Pennington 

Introduction.* 

If  two  different  fluids  having  a  common  plane  boundary  are 
accelerated  In  a  direction  perpendicular  to  the  boundary,  any  small 
Irregularities  In  the  boundary  will  tend  to  change  In  shape.  If 
the  acceleration  Is  directed  from  the  more  dense  to  the  less  dense 
medium,  the  Irregularities  will  tend  to  smooth  out  (in  the  absence 
of  external  forces).  Thus  the  plane  configuration  of  the  inter¬ 
face  is  a  stable  one.  This  can  be  illustrated  by  the  usual  example 
of  a  glass  of  water  sitting  at  rest.  If  one  considers  the  force 
of  gravity  to  be  replaced  by  an  acceleration  which  produces  the 
same  effect,  the  vater  and  the  air  are  undergoing  an  upward  acce¬ 
leration.  Since  the  acceleration  is  from  the  more  dense  to  the 
less  dense  medium,  the  air-water  interface  Is  stable. 

Returning  to  the  general  case.  If  the  acceleration  is  directed 
from  the  less  dense  to  the  more  dense  medium.  Irregularities  of 
the  Interface  will  tend  to  grow.  This  is  the  effect  known  as  Taylor 
Instability.  An  example  is  the  case  of  glass  of  water  turned 
upside  down.  Here  again  the  force  of  gravity  may  be  considered  to 
be  replaced  by  an  upward  acceleration.  The  acceleration  is  from 
the  air  to  the  water,  and  the  air— water  interface  is  instable. 


* 


The  results  In  this  paper  were  obtained  1..  1951  while  both  authors 
were  In  residence  at  Princeton  University. 
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The  water,  instead  of  maintaining  a  nearly  plane  lower  surface 
as  it  falls,  will  tend  to  Jet  out  into  long  spikes.  It  is  the 
formation  and  rate  of  growth  of  these  spikes  which  we  wish  to 
Investigate,  taking  into  account  the  effects  of  surface  tension 
and  viscosity. 

Section  1:  Taylor’s  Results. 

Let  us  begin  by  presenting  an  account  of  the  work  done  oy 
Taylor  himself  [4J  (see  also  [2]).  The  model  used  is  that  of  two 
fluids  of  infinite  depth.  The  Interface  (neglecting  perturbations) 
Is  the  plane  y  *  0,  the  y  axis  being  vertical.  The  initial  pertur¬ 
bation  will  be  of  the  form  cos  kx ,  with  amplitude  small  compared  to 
wave  length.  The  problem  is  then  two-dimensional,  and  the  true 
equation  of  the  Interface  at  any  time  is  y  »^(x,t),  where  the  func¬ 
tion  ")^(x,t)  is  to  be  determined  from  hydrodynamic  considerations. 

The  fluids  will  be  considered  to  be  incompressible,  and  only  the 
linear  terms  in  the  equations  of  hydrodynamics  will  be  used. 

The  linearised  hydrodynamic al  equations  in  either  fluid  are: 

(1)  *  vy  -  0 

(2)  ut  +  px  "  0 

(5)  »,  ♦  y  Py  +  B  +  Ji  ■  0 

Here,  as  usual,  u  -  component  of  velocity  in  the  x  direction, 
v  -  component  of  velocity  in  the  y  direction,  p  -  pressure,  - 
density,  g  -  acceleration  of  gravity,  and  gi  is  the  upward  accelera¬ 
tion  of  tue  system.  These  equations  have  solutions  of  the  form 
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(4)  u  -  -  4)x,  v  -  -  <t)y 

(5)  p  -  pD  -  (g+gJ/’y  +  '4>t* 

where  (J)  +  $)  -  0  and  pQ  is  the  mean  pressure  at  the  interface 

3cx  yy 

in  the  unperturbed  condition. 

For  the  upper  fluid  we  take 

(6)  -  Ae-*^  f(t)  cos  Kx 

(V)  Pi  -  PQ  -  (g+gi)/3^  +  A  ($i  )t 

tfnd  for  the  lower  fluid, 

(8)  $8  -  -  AeKy  f(t)  cos  Kx 

(9)  Pa  -  P0  -  (g+gi)/°2y  +  Z3 a  ( )  t 

the  lower  fluid  being  the  more  dense,  l.e.,^2  >  f*i '  ^he  al)0ve 
relations  satisfy  the  conditions  that  velocities  are  finite  at 
y  -  oo  and  y  -  -  oo ,  and  that  vt  -  v2  at  the  (approximate)  Interface. 
The  free  boundary  condition  at  y  -  7^(x,t)  is  that 

(10)  -jft-(y  “  >?(*»*))  -  0  or» 

<ll>  m  0  or' 

(12)  v  -  y[xu  -  /jt  ■  0 

neglecting  the  non-linear  term, 

(13)  y|t»v-KAf(t)  cos  Kx  or, 

t 

-  KA  f  f(t)dt  c os  Kx 


(14) 
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The  pressures  at  the  Interface  must  satisfy  the  condition 

(15)  Pi  »  p2. 

Substituting  from  (7)  and  (9),  we  have 

(16)  -  (g+gi  )/i  T]  +  fx  (4>i  )t  -  -  (g+g,  )/J2  7)  +/52(02)t  or, 

t 

(17)  -  (g+gt)(/2  -  A )AK  S  f(t)dt  cos  Kx  +  P2  (-Af'(t)cos  Kx) 

fco 

-  fx  (a  f'(t)  cos  Kx)  -  0  or, 

(18)  -  (g+gi  )  (  f>z-  /i  )AK  }  f  ( t  )dt  -  (  pz  +  fx  )A  f ' (t)  -  0 

fco 

or,  differentiating  with  respect  to  t, 

(15)  -  (g+gi)(/a  -fx)  Kf(t)  -  (fa  +  /?1)f"(t)  «  0 

so  that  we  may  take 

(20)  f(t)  -  sin  h  nt 

(this  choice  makes  the  fluid  velocity  zero  at  t  «  0).  From  (19), 

(21)  n2  -  -  _(g+gi)  (  P*  ~  A  )K 

P*  +  ft 

and  the  interface  is  given  by 

(22)  7j  -  KAn  *  cos  h  nt  cos  Kx 

If  (g+gi)  Is  negative,  (21)  has  a  positive  root.  The  dis¬ 
turbance  grows  like  cos  h  nt,  so  the  motion  of  the  interface  is 
Instable.  This  Instability  exists  for  all  positive  k,  i.e.,  for 
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all  wave  lengths  of  the  Initial  disturbance.  Note  that  the  smaller 
the  wave  length  (■  -^— ) ,  the  more  rapid  the  growth  of  the  dis¬ 
turbance.  This  limits  the  use  of  the  above  result  for  arbitrary 
disturbances.  Since  the  differential  equations  used  were  linear, 
one  would  expect  to  discuss  an  arbitrary  disturbance  by  Fourier 
analysis.  Let  the  surface  at  time  t  -  0  have  the  equation 


(23) 


y  -  f(x) 


cos  Kx 


Then  at  time  t,  we  have 


CD 


(24) 


y  - 


cosh  nt  cos  Kx 


where  n  Is  essentially  \Jk  .  For  t  +  0,  the  series  (20)  will 

oo 

diverge  unless  the  convergence  of  ^  aj^  cos  Kx  Is  extremely  rapid, 
since  cos  h  grows  so  rapidly. 


flection  2:  Viscosity. 

The  effects  of  viscosity  on  the  arguments  of  Section  1  are 
clear  Intuitively.  Viscosity  is  not  to  be  expected  to  remove  the 
instability,  but  only  to  reduce  the  rate  of  growth  of  the  amplitude 
of  the  disturbance  for  any  particular  frequency.  The  amount  of 
this  reduction  fcr  small  wave  lengths  Is  rather  startling,  however. 

In  particular,  as  the  frequency  - >oo,  the  rate  of  growth  of 

amplitude  — >  0. 

The  model  to  be  used  here  Is  that  of  Section  1.  The  (linearized) 
equations  governing  the  motion  of  an  Incompressible,  viscous  fluid 
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are 

(1)  ux  +  vy  “  0 

(2)  ut  +  -jr  Px  "  ^fT  (uxx  +  uyy) 

(3)  vt  +  -^5-  Py  +  g  +  Si  -  -j-  (v^  +  v^) 

where  Is  the  coefficient  of  viscosity. 

These  equations  are  satisfied  by 

(4)  u  -  -  -  yr  v  -  -  $y  + 

(5)  p  *  P„  -  (e+8i )  Py  *  /W t 
provided  that 

<6>  »xx  +  Vy  ‘  °-  f  -  (^XX  +  fyy>  *  ft,  of.  QJ  ,  [2]. 

For  the  upper  fluid  we  take 

(7)  «  Ae-Ky'fnt  cos  Kx 

-mi y+nt 

(8)  Yi  =  Be  sin  Kx 

(9)  Pi  -  P0  -  (g+gi)/°iy  +  PA$ i)t 
wher*2 

(10)  mf  -  Kz  +  — 

*x 

For  the  lower  fluid. 
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(10 

g 

-  COS  KX 

(12) 

r. 2 

r  m*y4nt  sin  Kx 
-  I  <e 

(n) 

P2 

-  P0  -  +P2(l)2)l 

wne  re 

do 

-  KJ  +  ^  " 

In  order  that  the  velocity  components  remain  finite  for  y  — > 
It  is  necessary  tliat  the  real  parts  of  m ,  and  n\a  be  positive. 
Let  the  Interface  be  giver  by  y  -  7j(x,t).  Th.-n 

(15)  >Jt  -  Vi 

or 

( 1 6 )  /Jt  -  K(A+B)ent  cos  Kx 
from  which 

(17)  -  K(A4B)n"1ent'  cos  Kx 

The  boundary  conditions  at  tne  Interfax  are 
(10;  t>j  -  u2 

(  1  9 )  v  I  ■  V  2 

(20;  -  p,  +  2/<,  -|yi  -  -  p*  ♦  2/^ 

(21)  M,(  4tL  +  -iHl,  .  yua  (  ♦  -1^2.) 

v  /  dx  dy  /  dy 


+ 

-  OD  , 
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7’he  last  two  equations  3tate  the  equality  of  the  components  of 
the  stress— tensor. 

Substl  ti.tion  in  equations  (l8)  to  (2!)  gives  four  conditions 
on  the  constants  A,  B,  C,  and  P. 

From  ( l6 ) , 

(22)  KAent  sin  Kx  +  mtBe  1  sin  Kx  -  KCent  sin  Kx  -  m2Dent  sin  Kx 
or 

(23)  KA  ■+■  rn i B  —  l) C  4"  rn2P  ■  0  1 

From  (19), 

(  ?4)  KAent  cos  Kx  +  3ent  cos  Kx  -  -  KCent  cos  Kx  +  KDent  cos  Kx 


or, 

(2)  A  +  B+  C-  D-  0 


Using  (20)  we  have,  after  some 

(26)  j  _-(g+sx )  )K  _  />in 


8lnpllf 1  cation , 


-  2^K2 


A  + 


-Ig+g.iIt/Wt  )k  _ 


+  n  +  2/feK'JC  -  -  0 

Henceforth,  we  will  let  -(g+g,  )  (/^-fl  )K  -  &  .  Similarly,  using 
(21)  and  soar*  algebra,  we  obtain 

(27)  2/^^*  +  A(K*  +  m?)B  +  2/feKaC  -/<*  (K2  +  n|  jD  -  0 

Equations  (23),  (25),  (26),  and  (27)  are  linear  ana  homogeneous  in 
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A,  B,  and  0.  They  have  non-trlvlal  solutions  If  and  only  If  the 
determinant,  of  thp  coefflclenta  vanishes, 


(28) 


K 

%  -  ft 


m 


A (K2+m? ) 
£  -  Kmj 


2>K2 

P 2  n  +  2/4^  K 2 


fP2 


7^2  ( K*  tmf  ' 

—2^2  K-ITlp 


J 


This  equation  reducea  to 


(29)  +  f*3  jj'/*'  K  ♦  )  +  (/*» ^  +  mi  ) 

4-  4nK  [/*t  K  +/l2m2j|^BK  + 


Equation  (20),  together  with  (10)  and  (l4),  gives  a  polynomial 
equation  of  „pnth  'degree  In  n.  Since  the  roots  cannot  ne  directly 
determined.  It  will  be  more  profitable  to  avoid  rationalization 
and  see  what  'nformatlon  can  be  obtained  by  other  means. 

In  Section  1  we  found  that  n  was  positive,  which  Implied  that 
Instability  occurred,  when  (g+gj  )  was  negative.  There 

(30)  n2  -  -tgy tVp -?>)*■_ 

p  +  r2 

The  value  of  n  which  determined  chiefly  now  fast  the  amplitude 
of  the  disturbance  grew  was 


(31) 


n 


-(g +zt)(pi-pxjy- 

f  I  f  2 
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In  th**  present  case,  then  expect  that  when  t-  -(g+gi  )  (Pz-fi  )K 
Is  p  jp  1 1  I  vp  ,  then'1  will  be  at  lead  jtp  r  ot  of  equation  (29)  with 
positive  real  part,  and  that  for  ih  Is  root, 


(32) 


Re(n)  £ 


/  S 

V  T/Pi  4  /’a  ) 


/-(g+g i  )  {Ps-  P\  )K 

/i  4/>a 


We  shall  determine  whether  equation  (2^ )  has  a  root  with  posi¬ 
tive  real  part  by  considering  n  a  complex  variable  and  applying 
the  principle  of  the  argument  to  the  right  half  plane.  So  long 

an  n  remains  In  the  right  half  plane,  the  quantities  rr.j  ■  \J)\3  + 

/ - T~n —  ^ 

and  mo  »\/K2+  ~ —  remain  on  one  branch  of  their  domain  of 

v  ^ s 

values,  so  we  will  run  Into  no  confusion  If  we  write  equation  (29) 
as 


(35)  r.0  +  j/«,K  +^<g  n  +^2K  *  Ml  fi  nj 


;'nK  /SK  fi n  K*  +/*,/°,n  - 


0 


Let  us  use  the  contour  C  consisting  of  the  part  of  the  Imaginary 
axis  between  (0,R)  and  ( 0 ,  — H )  and  the  Bemi-clrcle  in  the  right  half 
plane  with  this  segment  as  diameter. 

If  we  denote  t.ne  left-hand  side  of  equation  (3;)  by  f(n),  the 
principle  of  the  argument  states  that 


(34) 


1  [change  In  argument  of  f(n)  around  (T] 
2tt1  [number  of  zeros  of  f(n)  within  C] 
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ao  that 


(35) 


2Tl£number  of  zeros  of  f(n)  In 
right-half  plane] 


Now 

(36) 


f ' (Re 


10, 


^(He1^ 


1  Rp 


+ 


f'(n) 

nnp 


i:. 


provided  that  the  limits  on  the  right-hand  side  exist.  Since  they 

do,  we  shall  evaluate  them  separately.  To  evaluate  the  first  term, 

h  /2 

note  that  the  highest  power  of  n  appearing  Is  n  '  .  In  the  limit 

this  Is  tbs  only  term  which  will  matter,  so 


(37) 


llm 

R - >  oo 


ReUd0 


_  /r\  n 


-  llm 


oo 


ir/2 

s 

—w/2 


Re‘V/2 


(Re 


^7572" 


!Re1J)d0 


v/2 


S 


5/2  ld0  «  5/2  lir 


/2 


Hence,  In  the  limit  the  change  of  argument  of  f(n)  when  n 
traverses  the  semicircle  Is  5u'/2.  The  second  term  In  equation  (36) 
la  1  times  the  change  of  argument  of  f(n)  as  n  traverses  the  Imagi¬ 
nary  axis  from  +1  00  to  -1  00 .  This  change  of  argument  can  be  seen 
directly. 
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Comblning  the  various  changes.  It  Is  not  difficult  to  see 
that  change  In  argument  of  f(n)  as  n  traverses  the  Imaginary  axis 
is  from  5v/J»  to  3tt/4  or  -tt/2. 

Hence 


(38, 


5tt1/2  +  ( —1  /2t  )  1  -  2rrl 


Therefore,  f(n)  has  one  zero  In  the  right  half-plane,  which 
means  that  instability  does  occur  for  /3  >  0. 

We  next  ask  for  an  upper  limit  on  the  value  of  the  real  part 
of  this  root.  We  note  first  of  all  that  the  root  itself  is  real. 
For  positive  reaj.  n,  f(n)  Is  real  and  continuous  in  n  for  all  k. 
For  n  -  C,  f(n)  -  —  2  +/*z)K  <  0.  For  n  ■  \Z$/f 1+^*2  » 


f(n) 


4nK  Qii  K  +V^|K2  +  N/£?K2  +^n]  >  0. 


Hence  there  Is  a  positive  real  root  between  n  -  0  and  n  ■  +  A 

We  already  have,  then,  an  upper  limit  on  the  root.  At  the  root, 
n  <  +  •  This  is  Just  as  we  expected,  for  the  positive  value 


of  n  when  viscosity  was  neglected  (Section  1,  Equation  (l4))  was 

We  can  find  an  upper  bound  on  this  root  which  shows  more  about 
the  nature  of  the  root.  To  do  this,  we  rewrite  Equation  (33)  as 


(39)  |  (g»  +g)  {fc—px  )K  +  (Pi+(°2  )n2l  j  - — -  1  - -  +  _  1 _ 

J  |_ A K+>§K2+^2yflfen  K+^ufK^/S /h n 


+  4nK  -  0 


P—403 

-13- 


(bear  In  mind  that  g  +  gt  la  negative).  Consider  alao  the  comparlaon 
equation 

(Oo)  [(g^,)(P2-p,)K+(f1  +  ft)zJl 

Por  any  K,  the  positive  root  z  of  (40)  must  be  greater  than  the 
positive  root,  n  of  (39)  •  The  second  factor  of  the  first  tenr  has 
been  increased,  aid  at  a  root  this  must  be  counter  balanced  by  an 
increase  in  the  second  terra  or  a  decrease  In  the  first  factor  of 
the  first  term,  both  of  which  require  an  Increase  in  the  root.  For 
an  upper  bound  on  n,  then  we  have  only  to  give  the  value  of  z.  We 
rewrite  (40)  as 

(*D  (fi  +pt)T?  +  2zK2  (/!+/*)  +  (g  +  gi)(f2  —  p.)  K  -  0 

The  positive  root  Is 

(42)  z  -  -(M»  +Ag)Kg  +  >/(^t  -*-^2)aK4  -(g+giMpg-ftHfa+fijK 

A  +  f 1 

The  most  Interesting  thing  about  this  root  Is  the  fact  that 
It  has  a  maximum  for  some  K.  Thus  the  introduction  of  viscosity  has 
eliminated  the  tendency  for  disturbances  of  small  wave-length  to 
Increase  without  bound,  ke  would  like  to  know  the  value  of  z  at  ths 
this  maximum,  since  this  value  will  be  an  upper  bound  on  n  for  all  K. 
Differentiating  (42)  with  respect  to  K  and  setting  dz/dK  =  0,  we 
obtain 


r _ , _  t 

/*>k+  k*+o  j 


i+  4zK 
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(43)  K  -  --U+gi)  (Pa-ft) 

4z^!  +yU2) 

Substituting  this  '  alue  in  (4l),  the  result  is 

(44)  z  .  £-(g+gi>  fo-P.]271 

This  occurs  for 

(45)  K  -  £lIg+gl)(/W» ifl ^l/3(/pi-*-/°g  ) 1/3 

so  that  for  all  K  we  have 

(46)  n  <  £-t£+8i  h/w, 

2(/l*/fe)1/7(/*,4i^)T7;r 

'One  can,  of  course,  make  better  approximations  for  n.  For 
the  general  case  this  process  does  not  seem  to  offer  much,  since 
the  general  state  of  affairs  is  now  established. 

A  quite  complicated,  but  straightforward ,  calculation  shows 
that  n  has  only  one  maximum. 

Note  that  in  this  present  Section  one  cannot  satisfy  the  con¬ 
dition  that  tne  velocities  be  zero  when  t  -  0.  Apparently  because 
of  the  linearization  performed,  one  obtains  no  motion  at  all  if 
one  attempts  to  satisfy  this  condition. 
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Sectlon  3-  Surface  Tension. 

We  now  Introduce  the  effects  of  surface  tension  Into  the 
arguments  of  Section  1,  [2].  It  Is  to  be  expected  that  the  pres— 

•nca  of  surface  tension  will  remove  the  Instability  for  suffi¬ 
ciently  small  wave  lengths.  This  Is  Indeed  the  case,  as  will  be 
shown . 

To  Introduce  surface  tension  Into  the  arguments  of  Section  1, 
we  merely  replace  equation  (15)  of  that  Section  by  the  condition 

(1)  />.  -  />,  T,  -  0 

Substituting  from  equations  (7),  (Q),  and  (21)  of  Section  1, 

we  have 

(2)  —  (g+gi  )  (^2— fi  )  +  /^(O  *)t  —  /?i(^i)t  +  Ti/^XI:«0  or, 

(3)  -  ( g+g i )  ( A -ft  )  A  Kn-1  slnh  nt  cos  Kx  -  (/'2+/;i)An  slnh  nt  cos  Kx 

-  Tj  AK3  n~ 1  slnh  nt  cos  Kx  -  0 


so  that 

(4)  n2  -  )  (/jgr-fi  1 _  K _ li _  K3 

A  ♦  f 2  A  +  A* 

The  condition  given  for  Taylor  Instability  was  that  g  +  gi  be 
negative.  But  we  see  from  equation  (4)  that  the  amplitude  of  the 
Initial  disturbance  grows  only  when 


T i  K3 


>  0 


(5) 


-  (g+gi)(A-ft ) 

A  +  A 


K 


ft  ♦  f  2 
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or 

(6)  K  < 

or 

(7)  A  > 

2V 

where  A  -  —  la  the  wave-length  of  the  Initial  disturbance. 

Thus  for  wave-lengths  smaller  than  those  satisfying  condition  (7), 
there  is  no  Instability! 

Another  fact  of  Importance  is  expressed  by  equation  (4).  Since 
the  right-hand  side  has  an  absolute  maximum,  there  is  a  "most  dan¬ 
gerous  frequency,"  i.e.,  a  frequency  for  which  the  amplitude  of 
the  disturbance  grows  most  rapidly. 

The  most  dangerous  frequency  is  that  frequency  for  which  n,  or 
n2  ,  is  a  maximum.  At  this  frequency,  then. 


(9)  K  -  — - - ./-  (g+gi  )(f»-  fi  ) 

fyrT  V  ' 


• 

This  explains  the  hanging  of  water  droplets  on  the  underside  of  a 
horizontal  surface,  such  as  a  celling.  Such  a  droplet  is  under¬ 
going  an  upward  acceleration  of  980  cm/sec2  and  will  tend  to  drip 
because  of  Taylor  instability  unless  its  effective  wave-length  is 
too  small  tJ  satisfy  (7).  For  water,  the  critical  wave-length  is 
about  A  *  2tt  J 7 4/890  *»  1.73  cm.  Droplets  of  larger  diameter  will 

tend  to  drip,  while  smaller  ones  will  tend  to  hang.  (Actually, 
of  c curse,  the  true  critical  diameter  will  be  different  because  of 
circular  symmetry,  etc.,  but  the  above  at  least  contains  the  prin¬ 
ciple  Involved.) 
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Substitutlng  this  value  In  equation  (4),  we  have 


(10)  n2 _ L_  [-(^8. )  >]  3-2 

j  lyf\  p*  +  pi 


It  Is  remarkable  to  note  the  small  effect  which  the  numerical 
value  of  the  surface  tension  has  on  the  rate  of  growth  of  ampli¬ 
tude.  Although  It  Is  the  quantity  which  places  a  limit  on  the  rate 
of  growth  of  amplitude.  It  Is  felt  numerically  only  In  the  one- 
fourth  power,  as  equation  (10)  shows. 


Section  4:  Viscosity  and  Surface  Tension. 

In  this  Section  we  combine  the  results  of  the  two  preceding 
to  give  an  over— all  picture  Including  both  surface  tension  and 
viscosity.  We  would  expect  that  as  In  Section  3*  there  would  be 
no  Instability  for  small  wave  lengths;  and  that  for  longer  wave 
lengths,  the  rate  of  growth  of  amplitude  of  the  disturbance  will  be 
less  than  that  given  In  Section  J>. 

The  procedure  will  be  to  take  the  arguments  of  Section  2,  where 
viscosity  Is  considered,  and  alter  them  to  Include  the  effects  of 
surface  tension.  To  do  this  we  must  replace  equation  (20)  of 
Section  2  by 


(1) 


“P*  +  2/Xa  +  p‘ 


2/S 


0Vl 

ay 


o 


Substitution  In  this  equation  from  equations  (7)  -  (17)  of  Section  2 
yields 
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(2) 


[-  -H -  +  f-n  +  2^.K2]a  + 

+  f-  2yU2Ka  -  />*nj  C 
where  «  -  p(g+g-  )  (fa-fi  )K  -  TjK3 . 


+  2 2  Knia  D 


0, 


The  other  three  conditions  on  A,  B,  C,  and  D  are  the  same  as  those 
in  Section  2,  namely: 


(3)  A+B+C— D«0 

(equation  (25)  of  Section  2) 

(^)  KA  +  iDj  B  —  KC  +  maD  «■  0 

(equation  (:3)  of  Section  2),  and 

(5)  2/t1KzA  +  /S  (Ka  +  mf  )B  +  2/<-2K8C  -/<*(K8  +  mf)D  -  0 
(equation  (27)  of  Section  2) 


Equations  (6)  —  (9)  are  linear  and  homogeneous  In  A,  B,  C, 
and  D.  They  have  non-trlvlal  solutions  If  and  only  If  the  deter¬ 
minant  of  the  coefficients  vanishes. 


(6) 


1 

K 


2/4  Ks 


1 

/*i  (K8+m?  ) 


1 

-K 

2/*Ks 


“  *  +  fm+2/SK*  _  ?  +2^Kmt  -2/**K2-^2n 


1 

mg 

(k8+«4I  ) 
2fa  Kma 
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This  equation  reduces  to 

( 7  )  A  +  { Pi  +  Pz  \  +  /l2~*  )  +  (^X  +  y^itnx  ) 

*  Mr  K  |^XK  •►/(araa  |  j^ijK  +  /*x®i  j  -  0 

This  Is  precisely  equation  (29)  of  Section  2  wltn  (i  -  -(g+g i  )  )K 

replaced  by  a.  -  —  (g+gi  )  (^2—  P\  )K  —  TiK3.  In  Section  3,  where  sur¬ 
face  tension  alone  was  considered,  we  found  stability  for  oC  <  0 
and  Instability  for  a  >  0.  We  shall  show  that  these  conditions  still 
hold . 

For  >  0,  the  result  la  lamed  late  from  Section  2.  Para¬ 
phrasing  the  results  of  Section  2  for  cC  >  0  Instead  of  0,  we 

have:  for  <  >  0,  equation  (7)  (where  f*\  n  and 

-^gKa+/<*fsn  )  has  Just  one  root  with  positive  real  pa.-t.  This 
root  Is  Itself  real  and  is  less  than  +/3a  •  We  will  return 

to  the  problem  of  a  better  estimate  of  this  root  after  proving 
stability  for  oc  <  0. 

To  establish  stability  for  0(  <  0,  we  again  apply  tne  principle 
of  the  argument,  as  in  Section  2.  The  result  Is  established  by  a 
series  of  straightforward  but  laborious  arguments  which  we  shall 
omit. 

We  turn  now  to  the  Instable  case,  aC  >  0.  We  found  that  in  this 


case  the  equation 
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( 1 ' )  I  ♦  ( P\  +  P*  .  n2  I 


.!  -  -7 

+  />?**  ♦  /*  i  Pi  n  4 


/SK  +Jj%¥?  +  /*-2  f*n  /*sK  +Jp\K*  +  /Aif)i 
*-  4nK  =  0 


had  one  positive  root. 

There  are  two  Immediate  upper  bcjnds  for  this  root.  The  first, 
already  given.  Is 

(9)  -  r.r  *  -  -  T,KS 

V  A  +/°e  V  /°,  +  f>* 

and  for  all  K, 


Relations  (9)  and  (10)  state  merely  that  the  rate  of  growth 
when  both  viscosity  and  surface  tension  are  considered  Is  less  than 
Itoat  lfrien  surface  tension  alone  Is  considered. 

The  second  upper  bound  on  n  comes  from  comparison  of  (8)  with 
equation  (40)  of  Section  2.  Since 

(11)  a  -  -(g+gi)(/>2-/°i)K  —  TiK3  <  (g+gi)(f8-ft)K, 

the  root  of  (8)  must  be  less  than  that  of  equation  (39)  of  Section  2 
for  given  K.  This  may  be  seen  In  the  following  way.  Suppose  the 
value  of  £  in  (3)  is  increased.  The  first  factor  of  the  first  term 
tends  to  become  mo.  e  negative.  A’-  Increase  in  n  will  decrease  both 
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factors  of  the  first  term  and  Increase  the  second  term,  to 
counterbalance  the  change  In  flC.  Thus  the  root  of  equation  (40) 
of  Section  2  Is  an  upper  bound  on  the  root  of  (8).  This  Is  merely 
a  statement  of  the  physical  fact  that  the  rate  of  growth  when  both 
viscosity  and  surface  tension  are  considered  la  less  than  that 
when  viscosity  alone  le,  considered. 

From  the  study  made  of  equation  (39)  In  Section  2  we  can  give 
an  upper  bound  for  the  root  of  (8),  namely, 

(12)  n  <  ~(A  +  AgjK2.*  ...  +/*a^4  ~  (g+Sl  )  (f?-A  )  JEIlfl 

(f*  +  Pi  ) 


and  for  all  K, 

(13)  n  <  H  2/3 

*(Pi+hV/7(/H+fi)'7* 

The  upper  bounds  on  n  given  by  (9),  (10),  (12),  ar.d  (13)  will 
not  usually  be  of  great  practical  value.  For  particular  cases, 
numerical  methods  must  be  used. 

A  little  can  be  aald  about  the  frequency  for  which  (8)  has 
maximum  root.  The  effect  of  viscosity  Is  to  shift  the  maximum  toward 
smaller  K,  or  greater  wave  lengths.  Furthermore,  n  has  a  unique 
maximum  as  a  function  of  K. 
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Sectlon  Numerical  Examples. 

In  order  to  demonstrate  the  effects  of  surface  tension  and 
viscosity,  we  give  some  examples  for  ordinary  fluids. 

Example  1.  If  the  two  fluids  Involved  are  air  and  water,  sur¬ 
face  tension  would  be  expected  to  play  an  important  role  In  the 
development  of  Taylor  instability.  We  use 

(1 )  f air  -  0 

( 2  )  f water  -  1  g/cc 

(3;  ?i  -  74  dynes/cm 

(4)  g  ♦  gj  -  -2  x  104  cm/sec2  ^  -20 g 

Figure  l  shows  values  or.  n  vs.  k  when  surface  tension  is  con¬ 
sidered  and  when  it  1b  neglected.  The  corresponding  equations  are: 

(5)  na  -  2  •  104  K  -  74K3  wd, 

(6)  n*  -  2  x  I04  K 

For  the  surface  tension  case,  n  has  a  maximum  of  about  355 
at  K  -  9-5  (  P  -  0.66  cm)  and  drops  to  zero  at  K  -  l6.4  (  P-  O.38  cm). 
The  deviation  from  the  no-surfac"  tension  case  is  indistinguishable 
for  K  <  3  (  P  >2.1  cm)  . 

Experiments  have  been  made  by  Lewis  [3]  for  accelerations  on 
the  order  of  that  used  above,  at  wave-lengths  on  the  order  of  one 
centimeter.  However,  the  published  results  are  not  in  a  form  which 
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allow  comparison  with  those  given  above.  It  would  appear  that 
experiment* 1  verification  of  the  effects  of  surface  tension  should 
not  be  difficult  to  obtain  with  apparatus  like  that  used  by  Lewis. 

Example  2 .  If  the  two  fluids  Involved  are  air  and  glycerine, 
both  surface  tension  and  viscosity  would  be  expected  to  play  an 
Important  role  In  the  d  velopment  of  "'aylor  Instability.  We  use 

(7)  pair  =  0 

(8)  ^glycerine  *  1.26  g/cc 

( 9)  /*alr  *  0 

(10)  ^-glycerine  -  14.9  poises 

(11)  Ti  -  63  dynes/cm 

(12)  g  +  gi  ■  2*104  cm/sec2 

Figure  2  shows  values  of  n  vs.  k  under  four  different  condi¬ 
tions: 

1.  Neither  surface  tension  nor  viscosity  acting. 

2.  Viscosity  only  acting. 

3*  Surface  tension  only  acting. 

4.  Both  viscosity  and  surface  tension  acting. 

In  this  way  the  relative  Importance  of  the  two  effects  for 
various  wave  lengths  are  made  apparent.  The  corresponding  equations 
are: 


n 


,a  - 


(13) 


1. 


2  * 104K 
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(14) 


(15) 

(16) 


2.  P2*  104(1.26)K  +  1.26n2“!  j  ■- - 1  —  ...  +  — l — 

L-V/(l4*S)iKi?+(l4.9)  (1.26)n  14.9K 

+  4nK  =  0 

3.  n2  =  2 •  104K  -  K3 

4.  r.2- io4(l  .26 >K  +  63K3  +  1 . 26n*  |  - —  --  —  1  ■  -  ■ 

-\/(l4-9)2K2+(i4.9)(l.26)n 


1 4 . 9K 

+  4nK  *  0 


It  Is  seen  that  the  viscosity  is  unimportant  for  K  <  1  ( ?\  >  6.28  cm) 
and  that  the  surface  tension  Is  unimportant  for  K  <  3  (^  >  2.1  cm). 
Experiments  have  been  made  by  Lewis  [3]  for  accelerations  on  the 
order  of  that  used  above,  at  wave-lengths  on  the  order  of  one  centi¬ 
meter.  It  would  seem  that  the  viscosity  effects  would  be  apparent 
In  these  experiments.  Thio  would  lead  to  an  observed  value  of  n 
much  smaller  than  that  preducted  by  the  theory  for  non— viscous 
fluids.  However,  the  experiments  gave  an  observed  value  n  greater 
than  that  predicted  by  the  simple  theory.  Lewis  explains  this  on 
the  oasis  of  viscous  drag  on  the  channel  sides  In  the  apparatus. 
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